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weakly singular Fredholm integral 
equation in reproducing kernel spaces 
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Abstract 


In the present paper, we propose a method to solve a class of weakly 
singular Fredholm integral equations of the second kind in reproducing ker- 
nel spaces. The Taylor series of the unknown function is used to remove 
the singularity and bases of reproducing kernel spaces are used to solve this 
equation. Efficiency of the proposed method is investigated through various 
examples. 
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1 Introduction 


The Fredholm integral equation with weakly singular kernel arises in differ- 
ent problems of mathematical physics such as potential problem, variational 
equilibrium, fracture mechanics, infrared radiation, and elastic contact prob- 
lems [5,20,21]. Although in [01,116,722] authors obtained analytical solution 
for integral equations with weakly singular kernel in special cases, but gen- 
erally this is not an easy task. So the numerical analysis standpoint plays 
major role for solving such equations. 
In this paper, we consider Fredholm singular integral equation 
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po(x)u(x) + A(x) [. ee ay = f(z), |ja| <<land0O<a<1l, (1) 


in which p(x) 4 0 and A(x) ¥ 0, for all a € [—1,1], and y,A, f, and k are 
known smooth functions and wu is the solution of equation (M]) to be deter- 
mined. We assume that problem has a unique smooth solution u. Several au- 
thors considered the numerical solutions of Fredholm integral equations with 


weakly singular kernel. In [2] Jiang and Cui considered integral equation of 
a 
they solved the problem in reproducing kernel space W1(0, 1]. Chen and Zhou 
considered second kind Fredholm integral equation with Hilbert type singu- 
larity [RB]. They used transform to remove singularity and solved problem 
in W(0, 27] using reproducing kernel method. In [2] Babolian and Arzhang 
Hajikandi solved (0) with k(z,y) = 1. Du, Zhao, G., and Zhao, C. consid- 
ered integro-differential equation with logarithmic kernel and Kalman ker- 
nel with boundary values [0]. They used smooth transform to remove 
singularity, solving the converted equation with reproducing kernel method 
in W2[0, 1]. Chen and Cheng in [@] used piecewise homotopy perturbation 
method (PHPM), for solving integro-differential equation with weakly sin- 
gular kernel. In [8] authors solved (]), using Taylor series of the unknown 
function u to remove singularity, and then Taylor expansion of k& together 
with Legendre polynomials as bases to implement Galerkin method. The 
Sinc-collocation method is studied by Maleknejad, Mollapourasl, and Os- 
tadi, to solve nonlinear Fredholm integral equations with weakly singular 
kernel [15]. Beyrami, Lotfi, and Mahdiani solved Fredholm integral equa- 
tion of the second kind with Cauchy kernel [6]; they removed singularity by 
smooth transform and used reproducing kernel Hilbert space (RKHS) method 
to solve problem in W3/0, 1]. Nili and Dastjerdi in [8] solved weakly singular 
Volterra-Hammerstein integral equation with operational Tau method. We 
use Taylor series expansion of u(y), at point y = x, to remove singularity; 
then by use of reproducing kernel space W™[—1, 1] we convert this equation 
to a system of linear equations. We demonstrate the method with conver- 
gence rate O(h™); so when we increase m and move from one reproducing 
kernel space to other, the rate of convergence will increase. 


first or third kind with weakly singular kernel of the form k(x, y) = 


The rest of the paper is organized as follows. In section 2 we present 
definitions and some useful properties of reproducing kernel spaces. In section 
3 we remove singularity, by use of reproducing kernel space to implement our 
method. Section 4 is devoted to error estimate and convergence analysis 
of our method. Section 5 contains some numerical examples illustrating the 
application of the proposed method. We end the paper with some conclusions. 
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2 Preliminaries and notations 


In this section we briefly review reproducing kernel properties of the Hilbert 
space W™ a, b| and also fix notations used in this paper. The Hilbert function 
space, Wa, b], is defined as the linear function space 


W™ [a,b] = {fl f, f,..., fo" are absolutely continuous, f(™ € L?{a, b]}, 


which is equipped with the following inner product 


m-1 


b 
(fawm = >> fO@g%(a) +f PO) (a2) g™) (a2) dr. (2) 


4= 


The inner product (2) induces the following Hilbert norm 


IIfllwm = VF, fw. 


The following theorem presents an interesting property of the Hilbert space 
W'™ Ia, b]. 


Theorem 1. [9] The Hilbert function space W™|a,b] is a reproducing ker- 
nel space with the conjugate symmetric reproducing kernel R(x, y), that is 
given by 


IRm(x,y) 2 pal aly)x'}, a< Y 
Rm(2,y) = 
aey=yo6ye tly 


in which coefficients c;(y) and d;(y) are the solutions of the following system 
of differential equations 


gQ2™m xr 
(-1)"2 Rew) — §(@ —y), 


ua ‘ 2m—i—1 
Se ae (3) 
eS ONG G01 


where 0(x — y) is the Dirac delta function. 


Remark 1. [5] According to system (8), by use of the Dirac delta function 
properties, the equation 


ij a5 te 9) 


Ox2™ 


is convert to 
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Oo IRm(«,y) 
Ox? 


_ 2 Rm (2,) 


sae i=0,1,...,2m—2, 


xr=Yy xL=y 


(4) 


a"-1 IR, (a,y) 
Og2m—1 


a?" 1 rR (x,y) 
Ox2m™-1 


a=yt 
By solving the system (@) with boundary condition of (8) with Mathematica, 
the coefficients c;(y) and d;(y) are computed. 

In particular, each function f € W™[a, b] satisfies the following reproduc- 


ing property 


f(y) = (f(.), Rm(..y))wm Vy € [a, b]. 


Theorem 2. [S]Let R(x,y) be a reproducing kernel of W™[a, b]; then 


Os Riayy) 


2 : a -_ 
Du‘ Byi € L*a, }], t+j =2m-1 


with respect to x and y and 


at R(x, y) 


axioyh O0<i4+ 7 <2m—-2, 
x’ Oy 


are absolutely continuous functions in [a,b], with respect to x and y; so we 
have 


Ot) R(a, y) ae 
~Outoyk EW a, d], i+j=m 


Lemma 1. Let f be a smooth function of order m on [a,b]; then, for k = 
1,2,...,m, we have f € W* a, bj. 


Proof. Let f be a smooth function of order m on [a,b]; thus f‘”), for k = 
0,1,...,m, are continuous functions on [a,b]. Therefore f) € L?[a,b], for 
k =1,2,...,m, beside the real constants MM, exist which we have 


f(a) < Me Va € [a,8); 


trivially f, for k = 0,1,...,m-—1, are absolutely continuous functions, and 
proof is complete. 
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3 Removing the singularity and implementation of our 
method 


Throughout this section, we use the Taylor series to remove singularity of ({); 
then by use of reproducing kernel space property, we construct our method. 
With Taylor series of the unknown function u(y), at x € [a,b], we have 


—@ 2 ” —-f k 
u(y) = ula) + (y — aye (x) + Py" (a) 4-4 PE (a) 
—7 k+1 
, ea Goa (5) 


in which €,., is between x and y. Using (Q), thus equation (M) converts to 


k (2) | 
aute) +r) ( OI fo eeniy— 2) ey 
i=0 =e 
L (b+) 
tf ae ee. ny -oedy) = £0) 


and singularity is removed; now we are in position to define operator L : 
w™(-1,1]) > wW™[-1, 1, 


ul (.) 


i) 


k 1 
L(ul.)) = w(.)u(.) waa(S / KC. y)(y — .)* dy 
i=0 a 


1 4,(k+1) 
, [, ek y)(y = ye-eay) ; 


By use of operator L, the equation () is convert to Lu(x) = f(x). It is 
easy to prove that, for 0 < a < 1, operator D is a bounded linear operator on 
[—1, 1] and that, for a = 1, L is a bounded linear operator on [—1 + €, 1 — €| 
for « > 0. Let {x;}22, be a dense subset of interval [a, b], for i = 1,2,...; take 
y = %;, in R,(a,y), and define y;(.) = Rm(.,2;). Moreover, assume that 
W, = L*y;, where L* is the adjoint operator of L. 


Theorem 3.Let {x;}%, be a dense set in [— 


, 1]; then the functions wy; € 
W™(-1,1], and {y;}%, is complete in W™[-1, 1]. 


1 
1,1] 


Proof. By using the reproducing kernel properties of W™[—1, 1], we have 
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Y=Xi 


Therefore, we get 


2m—-1 ; 1 
1 ORn(.,8 a 
bi = Uti) Rm(., 24) + (wi) > il Oks) [ kense- 2) ds; 
i=0 s=ai Y — 


hence py € W™[-1, 1]. Let ue W™[-1, 1], andi € N. Then we have 


(u, Wi)wm = (u, L* pi)wm 


Since {x;}%2, is dense in [—1,1], we have Lu = 0. By the uniqueness of the 
solution of (Hl), we get u = 0. Hence, {%);}S2, is complete in W™[-1, 1]. 


i= 


As a consequence of Theorem 3, we can expand each function as the 
following. Let {x;}92, be dense in [—1,1], and let wu be the solution of (Ml). 
Since u € W™[-1, 1], we get 


Co 


u(x) = ye aw; (a). 


i=0 


In what follows, we use u, as approximate solution of equation (0) which is 
defined by 


Um,n (2) = S > ainbi(a), LE [—1, all 
i=0 
and by use of operator L, we get 
Sa; Lui(a) > f(v), (6) 
i=0 


and with inner product of equation (@), with y;, we have 


n 


yea (Ludi, ej)wm =f, ei)we- 


i=0 
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Thus by replacing ~ by =, the singular integral problem converts to the 
linear system of equations Ba = F’, in which the coefficient matrix B = 
[b;,;] = [Lyx (a;)], and the right hand side vector is F = [F;] = f(a;). 


4 Error estimate and convergence analysis 


In this section first the convergence of the method is investigated, and then 
the error estimate of proposed method is studied. 


Theorem 4. Let {;}9, be dense in [-1,1], and let {y;}S2, be the orthonor- 
mal basis of W'™|—1, 1] produced by Gram—Schmidt process on {y;,}2,. Then 
Um,n, the approximate solution of equation (0), and its derivatives are con- 
vergent uniformly to u; that is, the exact solution of equation (Ml), and its 
derivatives, respetively. 


Proof. With the Gram process we have 


k=1 


and hence we get 


ee ). 


Therefore the approximate solution of equation ({) is 


tina(nh= b> S- aiBinPr (x) 


i=0 k=1 


and clearly we have limps oo ||U — Um,n||w~[-1,1] = 0. On the other hand, we 
have 


m-1 


u(x) = (u(.), Rm(.,2)) = wOQRRG2) + f wmEROM, a)dz. 


i=0 


Differentiating 7 times with respect to x, we get 


m— 1 i p(m) 
- Sx 2 bis 0 x) +f ul) (2) OF Rin (2,2) 4 
=a Ox 
=0 
= OR Gt) 
= (u(y), her 


Thus we have 
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a a Otel ,&) 
ju (2) — wg) | = (a6) = mn em 
ORm(.,x 
<I] 2 hn wma] Mayen 


Since se € W™{-1, 1] and || OR) |lw~[-1,1], is continuous with 


respect to x, on [—1,1]. For some M > 0, we have 
OB (oy) 


mi_p ys M. 
| Fai lwrm-ays M, 


which completes the proof. 


Let {x;}"_, be a subset of [—1,1], and let -1 < a < a < ++: < 
Lyn < 1. We define hj := 2341 — @j, 7 = 1,2,...,n —1, and also put 
h:= Max1<j<n—-1 hj. 


Theorem 5. [AjLet u €¢ W™[-1,1] be a smooth solution, and let Um n be 
the approximate solution of the equation (0). Then 


||U — Um, nlloo < ch™, 


where c is a constant. 


5 Numerical examples 


We examine accuracy and convergence of the proposed method through four 
examples, which indicate efficiency of the method. The errors were defined 
as 


Emon = ||U- Um nlloo & max | u(t) — Unn(x:) | m=1,2,.... 
l<i<n 


Example 1. [8] Consider the following integral equation 


ap lees xy + 2y?)uly 
ula) + f | JMO ay = Fe) 
-1 (y— a)» 
where f is chosen such that 
u(x) = 2°, x € [-1,]], 


is the exact solution. The numerical results are given in Table I. In Table [I], 
for fix n, by increasing m, the absolute errors converge to zero. In Figures [I], 
for m = 6,8,10, when n = 8 the absolute errors of u and its derivatives are 
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Table 1: Values of Em,n, for Example 1. 


m n=8 n= 16 n = 32 

6 1.72675 x 107! | 1.34164 x 107! | 1.06141 x 107+ 
8 9.37476 x 10-4 | 3.09832 x 10-4 | 2.91592 x 1075 
10 || 6.73973 x 10-® | 2.6006 x 10-® | 3.41363 x 107° 

Table 2: Values of Em,n, for Example 2. 

m n=8 n= 16 n = 32 

4 2.21811 x 107? | 2.89879 x 10-7 | 7.6194 x 107? 
6 4.92002 x 10-4 | 7.14302 x 10-4 | 1.95649 x 107-3 
8 2.02376 x 107° | 1.07565 x 107° | 2.91592 x 107° 
10 |} 3.13506 x 10~-® | 1.2141 x 10-7 | 2.02656 x 10-7 


shown. 


Example 2. Consider the following integral equation 


1 aty ( 
e* uly) 
u(x) + o+2) | ————dy = f(a), 
(a) + (+2) f Tay = F(a) 
where f is chosen such that 
u(x) = ze”, x € [1,1], 


is the exact solution. The numerical results are given in Table 2. In Table 2, 
for fixed n, by increasing m, the absolute errors converge to zero. In Figures 
Q, for m = 6,8, 10, when n = 8 the approximate solutions and its derivatives 
have high accuracies. 


Example 3. [1,8,(13,14)Consider the following integral equation 


— ft uly) 
ae) rer 


where f is chosen such that 


dy = f(z), 


u(x) = @, x € [0,1] 

is the exact solution. For m = 4,6,8, when n = 8, the absolute errors of u 
and its derivatives are shown in Figures B. In Table B, the numerical results 
of proposed method are compared with [1,8 [13]. where the Max Error is 
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Figure 1: Example 1, m = 6(blue dotted), m = 8(red dash), m = 10(Gray line). 
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Figure 2: Example 2, m = 6(blue dotted), m = 8(red dash), m = 10(Gray line). 
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Table 3: The Max Error result of Example 3. 


n n= 16 n = 32 
Babolian’s method [3] 1.59 x 10° 
proposed method for m = 10 5.29781 x 107! | 4.75211 x 10~" 
Lakestani’s method [13] 127% 10-* 2713 10-* 
Product integration method [II] Ls t=" 9.39 x 10-7 
Lagrangian interpolant [I] 2.12 10-* 1.94 x 10~® 


Table 4: The Max Error result of u(.), Example 4. 


ao eo FS 


3.17985 x 107-4 
2.38288 x 10° 
2.93877 x 10~° 


4.44036 x 1074 
2.63628 x 1078 
7.58042 x 1079 


4.24008 x 10-4 
2.23585 x 10° 
3.42293 x 1078 


defined as Max Error = maxo<z<1 | u(%) — Umn() | - 


Example 4. [(3,19]|Consider the following integral equation 


1 


a) 
ue) — 35 ff Se t= FO) 


where f is chosen such that 


u(a) = 27?(1— <7), x € [0,1], 

is the exact solution. For m = 4,6, 8, when n = 8, the absolute errors of u and 
its derivatives are shown in Figures M. In Tables @, GB, and @ the numerical 
results of proposed method are given. where the Max Error is defined as 


Max Error = maxg<e<i | uw (x) — un (2) | . 


Table 5: The Max Error result of u“)(.), Example 4. 


com KTS 


2.10369 x 107+ 
8.89262 x 10~° 
3.66355 x 1077 


7.59437 x 1073 
2.43063 x 1078 
2.21809 x 1077 


7.61556 x 107? 
1.37019 x 10~° 
1.20238 x 10-7 
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Figure 3: Example 3, m = 4(green line), m = 6(blue dotted), m = 8(red dash). 


Table 6: The Max Error result of u‘?)(.), Example 4. 


n=8 


n= 16 


n = 32 


4.6866 x 1071 
1.18157 x 10-9 
1.3664 x 107° 


co oO FS 


3.45622 x 107! 
1.25931 x 1074 
2.79569 x 1077 


3.75916 x 1071 
1.20977 x 10-3 
1.67734 10-5 
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Figure 4: Example 4, m = 4(green line), 


(c) Absolute error for u?)(.). 


(d) Absolute error for u)(.). 


m = 6(blue dotted), m = 8(red dash). 
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6 Conclusion 


Different problems have been solved by researchers using reproducing kernel 
Hilbert spaces (RKHS). They assume the unknown solution belongs to spe- 
cial fix W3” and use Gram—Schmidt process to implement RKHS method. 
In this paper, we attempted to solve Fredholm integral equations of the sec- 
ond kind with weakly singular kernel. We used the Taylor series to remove 
the singularity. According to the problem, we investigated the problem in 
reproducing kernel spaces W™{—1, 1]. In the proposed method by use of re- 
producing kernel space property, the problem was converted to a system of 
linear equations. In the proposed method we do not use Gram process, we 
use RKHS property and use W5” for different m. Though we suppose the 
unknown solution is smooth, but this is a very strong assumption, and it 
is enough for the solution to be in W3” for some m. In this method for fix 
m, when the number of bases n increases the coefficient matrix of the linear 
system will be ill-conditioned; so we increase m and move from one reproduc- 
ing kernel to other instead, to find approximate solution with high accuracy. 
Several test examples were used to show the efficiency and applicability of 
the method. 
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